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In one dimension and for discrete uncorrelated random potentials, such as tight binding models, 
all states are localized for any disorder strength. This is in contrast to continuous random potentials, 
where we show here that regardless of the strength of the random potential, we have delocalization 
in the limit where the roughness length goes to zero. This result was obtained by deriving an 
expression for the localization length valid for all disorder strengths. We solved a non-linear wave 
equation, whose average over disorder yields the localization properties of the desired linear wave 
equation. Our results, not only explain the origin of the difficulty to observe localization in certain 
physical systems, but also show that maximum localization occurs when the roughness length is 
comparable to the wavelength, which is relevant to many experiments in a random medium. 


For more than half a century, thanks to the pioneering 
work by Anderson we have taken for granted that 
strong disorder will lead to localization of all states, par¬ 
ticularly in low dimensions and for uncorrelated disorder. 
Spatial localization is when an electron, an atom or even 
a photon, cannot propagate in a medium when disorder 
is large. Localization can also occur in time due to fluctu¬ 
ations [2] and was found to be relevant to the expansion 
of our universe [3] . Formally, the propagation probability 
decays exponentially with the medium’s length, which is 
known as Anderson localization. Fven in the presence 
of interactions between particles, strong localization is 
expected to occur, which is now popularized with the 
term many-body localization [4] . Most of the theoretical 
work, has focused on discrete random equations, such 
as tight binding models, where the theoretical results for 
Anderson localization are confirmed by numerous numer¬ 
ical studies and the main message can be summarized as 
’’all states are localized for uncorrelated disorder in one 
and two dimensions; in higher dimensions this is true for 
sufficiently strong disorder” 

Experimentally, localization has been observed in 
many different systems, including electrons [6], photons 
[7] and atoms [8l|9]. Moreover, applications are becom¬ 
ing increasingly important, particularly in random las¬ 
ing [ini E] and optics [T2]. These different systems all 
share a similar underlying wave equation. Here, instead 
of looking at discrete equations, we look directly at the 
continuous wave equation and show that for arbitrarily 
strong disorder, we have no localization when the rough¬ 
ness length vanishes, even in one dimension. While at in¬ 
termediate roughness, localization is maximized. This is 
in stark contrast to discrete models, where no equivalent 
delocalization occurs. To obtain this result, we used a 
new approach based on solving an equivalent non-linear 
disordered wave equation. Our result explains why it 
is sometimes difficult to observe localization in certain 
physical systems when the roughness length is not of the 
same order as the wavelength. 

Anderson localization {AL) has become an important 


phenomenon well beyond its original work on tight bind¬ 
ing models with random potentials and couplings, which 
describe quantum particles or spins [1]. AL is important 
in photonic systems HME] , random lasers [iniiii] , quan¬ 
tum information noise and entanglement |16j , atomic sys¬ 
tems 00, mechanical systems Bin], biological systems 
m, cavity QED m, as well as cosmology, where infla¬ 
tion is dependent on fluctuations [3]. All these systems 
share a common underlying wave equation, which can be 
written as 

[dl+p^{x)]i;{x) =0, ( 1 ) 

where tp is the amplitude, p{x) = -y/e — v{x), the clas¬ 
sical momentum, v(x) the random medium and e the 
energy. We will restrict our attention here to the quasi- 
one dimensional situation, where the effect of disorder 
is the strongest. However, many of these results can be 
extended to higher dimensions and will be discussed else¬ 
where. Solutions to equation 0 can be obtained for a 
random potential v(x) that is not continuous. For in¬ 
stance, if v(x) is written as a sum of delta functions or 
square wells, equation Q becomes equivalent to a tight 
binding equation studied by Anderson and others [Buni- 
[22] . The main result is the localization of all states if the 
potentials are uncorrelated, regardless of the strength of 
disorder. In the presence of correlations in the disorder, 
some states can be delocalized too |23U27| . It is impor¬ 
tant to note here, that the minimum correlation length in 
tight binding models is limited to the smallest distance 
between impurities or orbitals. However, when several 
next nearest neighbors coupling elements are non-zero or 
when there is mixing between different energy bands, this 
induces effective correlations between neighboring onsite 
potentials. Hence, in this case too, a continuous potential 
model is more adapted. 

When the potential v is continuous, the situation 
changes. For instance, we can consider a typical quasi- 
one dimensional (QID) random potential of width Ly. 
QID means that at low energies, we can restrict ourselves 
to the one dimensional wave equation 0 , where only the 
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FIG. 1. Top graph: a typical Gaussian disorder potential, 
v{x, y) with correlation length I connected by leads (in yel¬ 
low). For a quasi ID system in the lowest transverse mode 
the transverse potential can be integrated out to yield an 
effective ID potential, v(x) shown in the middle graph. The 
bottom graph shows the numerically calculated average trans¬ 
mission for such a random potential as a function of energy 
and correlation length 1. 

lowest transverse mode is relevant and the QID solution 
is simply sm(7ry/Ly)'ip(x), where ■i/' is the solution for po¬ 
tential v(x) = f sin(7ry/Ly)v(x,y)dy illustrated in figure 
We can consider the transport problem and evaluate 
numerically the transmission T, through such a poten¬ 
tial assuming that we have perfect leads or wave guides 
at each end and represented in yellow in the figure. The 
numerical result is obtained by discretizing equation 0 
and then computing the disorder averaged transmission 
for a given system length. Care is taken in choosing a dis¬ 
cretization parameter much smaller than both the disor¬ 
der correlation length, I, and the wavelength. This leads 
to the non-monotonic behaviour of the transmission as 
a function of energy and correlation length I shown in 
figure [2 At high enough energies and I, the transmission 
is maximum (1 in this model), while it is close to zero 
for a certain range of energies and 1. This is the strong 
localization regime (AL), which is usually discussed in 
ID random systems. In the opposite limit of vanishing 
correlation length I, the transmission is again maximum, 
which becomes a fully delocalized state at Z —>■ 0. We 
show below, that this regime is robust with increasing 
disorder as represented in figure 

For low disorder, we can understand the result in figure 
[l] using the perturbative approach to disordered poten¬ 
tials m- In this case it was found that when e 3> n, 
the inverse localization length, or Lyapounov exponent 
Am, is given by where Cy{x) = {v{0)v{x)) 

is the binary correlator of Fourier transform Cv{k), with 
ko = x/e the wavenumber, and (•) the disorder aver¬ 



FIG. 2. The Lyapounov exponent A as a function of corre¬ 
lation length I for different values of the disorder strength 
(here (j„ varies from 0.3 to 1.6 and e = 1.9, hence the poten¬ 
tial sometimes exceeds the energy). The binary correlator is 
taken to be Gaussian. The dots are the results obtained nu¬ 
merically for the decay of the transmission, with errors smaller 
than the size of the dots, while the lines are 2X^u from expres¬ 
sion (§ with no fitting parameters. The factor 2 comes from 
the difference in defining A from the transmission versus the 
wavefunction amplitude. 

age. This result leads to a delocalization-localization- 
delocalization dependence as a function of the disorder 
correlation length I shown in more detail in figure To 
compute A, we considered a Gaussian correlated poten¬ 
tial, i.e., {v{0)v{x)) = crle~^ with amplitude ay. 
Such a Gaussian binary correlator is obtained, for in¬ 
stance, when the potential, v{x) is a sum of Gaussian 
impurities located at random sites. Representative po¬ 
tentials with different correlation lengths I are shown in 
figure 0- The small disorder {ay e) result for the 
Lyapounov exponent is given by 

\ ^v{‘2'ko) _ 

^'^-8e~8e ’ ^ > 

which implies that for both I 0 and I —>■ oo. Am —>■ 0, 
while maximum localization occurs for I = lj2ko. In ad¬ 
dition, we also have Am —t 0 for e —>■ oo. This result 
differs substantially from the localization behaviour of 
the ID disordered Anderson tight binding model (AM), 
where localization occurs for all energies [5D]. Moreover, 
in the AM, localization (A) does not vanish at the small¬ 
est correlation length. However, long range correlations 
in the AM model can also lead to delocalization [551HH] j 
similarly to the continuous case shown here. In general, 
there is a decrease of the Lyapounov exponent with en¬ 
ergy as seen by the e~^ prefactor in equation 0, which 
is also true in higher dimensions [30] . However, correla¬ 
tions such as the roughness of the potential can override 
this behavior due to the exponential dependence on Z, 
which in some cases can even lead to delocalization at 
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small energies as illustrated in figure 
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FIG. 3. a) Examples of Gaussian disorder potentials with 
different correlation lengths {I — 0.01, I = 0.3, and I = 1) 
but with the same standard deviation ((j„ ~ 1.6). b) Corre¬ 
lation length (1) dependence of the Lyapounov exponent (A) 
for different values of the disorder strength {av from 0.07 to 
1.1) and e = 1.9. The dots represent A obtained numerically 
from the transmission, while the red curves (2Am) are from 
the perturbative expression given in equation ([^. 

At large disorder and when the correlation length is 
large, expression (1|) breaks down as seen in figures [| 
and |3l To understand the localization behaviour in this 
regime, which is relevant to many experiments, we need 
to go beyond the perturbative result, which brings us 
to our new approach to localization physics. The main 
idea is to solve an analogue to equation 0 but with an 
additional non-linear term: 

(3) 

Interestingly, there exists an exact solution to equation 
0, which can be expressed in terms of the integral so¬ 
lution 'ipix) = ip(0)e^-^o [31], where 

/(x) = p(a;) - [ K{x')e^^^’^^''>dx'. (4) 

Jo 
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Here Pi{x) = p{x')dx' is the integrated classical mo¬ 
mentum. If the right side term in equation (|^ vanishes, 
we recover our original equation 0 . Since this term is 
the difference between the classical and quantum momen¬ 
tum, we expect this term to be small and to vanish with 
disorder averaging (see appendix). Hence the localization 
behavior of the non-linear equation ([^ will describe the 
localization behavior of the linear equation Q. The last 
term /„ in equation 0 describes the memory effect of 
the wave propagation, expressed as an integral. For clar¬ 
ity, we have expressed the disorder dependence of p{x) 


in terms of ky{x) = p{x) — fcg with average 0 and vari¬ 
ance cr^ = {k1). fc(,(x) is the spatial derivative of kv{x), 
which we assume to be finite and which scales as l/l. We 
only consider the case where fc(,(x) remains finite, hence 
no discontinuous potentials. From here on, all the results 
will be expressed in terms of ky{x) rather than v(x). The 
reason is that at high disorder this is the relevant quan¬ 
tity, while at low disorder they are proportional, since 
v(x) ~ —2ky{x)y/t. The term fy(x) contains the physics 
relevant to the localization behaviour. Its average over 
all disorder configurations can be expressed in terms of 
a new correlation function Cp(x): 


{fv{x)) = [ (k'y(x')e~^"f^^' . 

Jo '' ---^ 


Cp(x,x') — Cp (x — 2 


( 5 ) 

The real part of the integral f{x')dx', which appears 
in the wavefunction solution, determines the wavenum¬ 
ber, while the imaginary part corresponds to the expo¬ 
nential dependence of the wavefunction. Hence we expect 
the disorder average of the imaginary part to be related 
to localization. Indeed, figure [^ shows the linear increase 
with x of fg '^{fy{x'))dx' with proportionality coefficient 
A. More precisely, we have 


-| pX poo 

A = S— / dx{fy{x)) 

^ Jo x^ao Jo 

( 6 ) 

assuming k^ real. This follows from the exponential de¬ 
pendence of the wavefunction, which determines A and 
can be expressed as A = ^ \n{\tf}{X)/'tlj{0)\)x^oo- Equa¬ 
tion (|^ is the main analytical result of this paper and 
is valid for all disorder strengths and correlations. It’s 
validity is illustrated in figure]^ For Cp{x) symmetric we 
have A = Qcp{2ko)/2. 

For arbitrary potentials, the correlator Cp{x) can be 
quite difficult to evaluate. However, it is possible to 
describe the localization behaviour in different impor¬ 
tant limits. The limit, where I is large and where the 
perturbative expression breaks down at large disorder, 
is important to understand for many experimental sys¬ 
tems, like in semiconductors in the presence of long range 
strong Coulomb potentials [35] or for scattering in pho¬ 
tonic crystals |33j . In this limit, where we consider a 
disorder potential characterized by a large but finite cor¬ 
relation length I and assuming that |x| <C I we can write 
(see appendix) 


Cpix) ~ ~ iPi2x)P;,{x^). 

( 7 ) 

We used that A:(,(0) is not correlated to any function 
of fc„(0), since (A:(,(0)fc"(0)) = 0 for any positive inte¬ 
ger n, when assuming random Gaussian impurities. The 
term Pp{x‘^) — {[{ky{x) — = 

—z(—ix)"“^([A:^(x) — fc„(0)]"+^) is related to the 
moment generating two point correlator and the Fourier 
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FIG. 4. Graph a) represents the Lyapounov exponent nor¬ 
malized by the disorder strength as a function of 1. The col¬ 
ored lines are the numerical results, while the black line is 
2Auj = c„(2fco)/4e from the perturbative expression (at the 
resonances Au, = 0). In b) the dots are from expression 0 
evaluated numerically, while the blue lines are the numeri¬ 
cal transmission results and the red lines are 2Au, from the 
perturbative expression, c) The distributions of the potential 
kv{x) (left) and k'^{x) for I ranging form 0.05 to 2. d) Plot 
of the increase of fgifv) with x (equation 4| shown in red, 
and the blue lines correspond to minus the wgarithm of the 
transmission —(InT(x)) evaluated numerically. A more de¬ 
tailed disorder strength dependence around the resonance is 
discussed in the appendix. In all figures, we used e = 40. 

transform of the distribution function of k'^ which sharp¬ 
ens for larger I as shown in figure |^. It can be written 
as Pp{x^) — 2[ck{x) — Cfc(O) -I- hm{x)\, where hm{x) are 
higher moments. In terms of notation, P'ik) = dkP{k). 
P{2x) = is the Fourier transform of the dis¬ 

tribution of ky, which is largely independent of I as seen 
in figure]^. The binary correlator, Ckix) = {kv{0)ky{x)) 
and consequently Cp{x) ~ 2i[ck{x) — Cfe(O)] in the lowest 
order of the disorder strength. This result can also be 
obtained directly from the Taylor expansion of the expo¬ 
nential term in Cp{x) and keeping only the first non-zero 
term. 

Assuming that k^l ^ 1, we can evaluate the Fourier 
transform of Cp{x) at 2fco, which gives rise to the following 
convolution: 

Cp(2ko) = 2i J dkP{k)[cl{2ko — 2k) + hm*{2ko — 2k)]. 

( 8 ) 

The first term is the distribution function of the disorder 
potential of width ak = \/{k^} centred at zero and nicely 
illustrates what happens with increasing disorder. For 
small enough disorder, P{k) is simply a delta function 
centred at fc = 0, hence using equation ([^ and dropping 
the higher moments, we have Cp(2fco) — 2jCfe(2fco) and 


A = Qcp{2ko)/2 ~ Cfc(2/co) — Ct,(2fco)/4e = 2Auj, which 
is twice the value of the perturbative result obtained in 
ref. [27] . This factor of 2 is due to the different aver¬ 
aging method. Indeed, in one dimensional disordered 
systems, we have for X —>■ 00 , A = ln(|'!/)(X)|)/X = 
2(ln |'0(X)|)/X = 2Au,, because |t/’(X)| follows a log¬ 
normal distribution [114] . Hence our result is equivalent 
to averaging the wavefunction directly. For larger disor¬ 
der the Lyapounov exponent becomes the convolution of 
the low disorder value at 2fco averaged around 2fco ± Scr^. 

This is best illustrated in the context of a disorder po¬ 
tential, where the Fourier transform of the binary corre¬ 
lator, has resonances. Such a potential can be obtained, 
for example, by starting with an uncorrelated potential 
and then making it smooth over a length scale I (see de¬ 
tails in the appendix and ref. |3S]). This leads to strong 
minima or resonances in the Fourier transform ampli¬ 
tude as reflected in figure]^ Precisely at the resonance 
{I = 0.85 in figure 4), the perturbative approach gives 
A„ = 0. With increasing disorder, the resonances, where 
A is minimum, first broaden (described by a convolu¬ 
tion) then reach a high disorder regime. At high disorder 
the full correlation function Cp(2fco) needs to be evalu¬ 
ated, which involves all moments. The second moment 
(or binary correlator) Cfe(2fco) |fc«(2A:o)P is simple to 
compute, since it is proportional to the squared absolute 
value of the Fourier transform of the disorder potential. 

In the other limit, where the correlation length I is 
small, we can see in figures and 0 that the behavior 
is largely independent on the disorder strength (except 
for the multiplicative factor). We can understand this 
result by first looking at the behavior of Cp{x) in the 
limit |a;| ^ , where 

Cp{x) ~ ~ iP^if)Pi2iy^\^K (9) 

For I —)■ 0 the exponential decay of Cp{x) vanishes and 
Cp{x) ~ iPp{P —>■ 0). For most disorder distributions we 
can write Pp{l^ —>■ 0) ^ —erf (for a Gaussian distribution 
the proportionality coefficient is one and P(2Z)I“I/^ —?> I 
for 1^0). Hence, to determine the Fourier transform 
of Cp{x) when Z —>■ 0, we can consider Cp{x) ~ —ierf for 
\x > I and Cp(a;) ~ —ia^x'^/P for \x\ < I using equation 
([S). This leads to Cp(2fco) — for Z —>• 0 (A ~ Zerf) 
and the result is valid for any disorder strength and only 
the proportionality coefficient will depend on the disorder 
distribution. Therefore, for arbitrary disorder strength, 
A will vanish linearly with vanishing Z. This delocaliza¬ 
tion can be understood, as the zero average of the disor¬ 
der potential within a wavelength. On the other hand, 
localization is the strongest when the wavelength is com¬ 
parable to Z and then A decays again at large Z. For very 
high disorder (a » e), we find numerically that A ^ a^l 
still applies for vanishing Z but then remains constant for 
larger Z k^^. 

To conclude, the (de)localization behavior has impor- 
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tant implications to our understanding of low dimen¬ 
sional systems. Often it is assumed that in one dimen¬ 
sion, all randomness will localize, but as we have shown 
here this is only the case for k^l ~ 1. For instance, in 
widely studied systems, such as GaAs bases heterostruc¬ 
tures, the disorder correlation length can be of the order 
of I ~ lOOnm, while the Fermi wave length is only about 
~ lOnm [36]. Hence we expect localization effects to 
be strongly suppressed when dominated by long range 
disorder. Depending on the disorder correlation, this 
suppression can be exponential (for a Gaussian binary 
correlator) or quadratic for an exponential correlator. In 
atomic systems this effect is important too, since usually 
I > lOOnm and the atomic de Broglie wavelength can 
be very small jSj. In the other extreme, of very short 
range disorder, like alloy scattering, the disorder corre¬ 
lation length is of the order of I ~ O.lnm (the atomic 
distance). For a typical Fermi wavelength of lOnm, this 
leads to an increase of the localization length by two or¬ 
ders of magnitude. A similar situation arises in photonic 
systems, where the wavelength is of the order of SOOnm, 
but if the disorder correlation length is much smaller, 
then no localization can be observed. We believe that 
the suppression of localization A ^ la'^ at small I is not 
necessarily unique to the continuous potentials we con¬ 
sidered here, but is likely to occur in other systems too. 
For instance, the equivalent tight binding model, with 
discretization a, would renormalize the disorder poten¬ 
tial by a^v for unit bandwith, which suggests delocaliza¬ 
tion for small a. More generally, any potential with fixed 
ay but vanishing integral over the wavelength isdikely to 
lead to delocalization. 

Summarizing, we have shown that the localization be¬ 
haviour of the standard disordered wave equation can 
be computed for all disorder strengths and correlation 
lengths using the disorder average of an approximate non¬ 
linear wave equation. This has important implications on 
our understanding of disordered systems and its applica¬ 
tions [37] as well as cosmological fluctuations [3EH]. 
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APPENDIX 

Derivation of equation (4) 

Equation (3) reads 

[dl+P^{x)]il’{x) = - {-id^ - p{x))'il){x)f 

( 10 ) 

taking ^ /(^ )<ta: leads to 

idxfix) + 2p{xf - 2p{x)f{x) = 0 (11) 

Assuming f{x) = p{x) + p(x) we have the following first 
order differential equation for 77 , 

dx[r]{x) + p{x)] + 2ir]{x)p{x) = 0, (12) 

with solution for 77 ( 0 ) = 0 and Pi{x) = p{x')dx', 

, 7 ( 0 ;) =r e^^^‘^^'Hd,,p{x'))dx'. ( 13 ) 

Jo 

Defining p{x) = ky{x) + fco we obtain equation (4). It 
is important to note that no approximations have been 
made beyond considering equation (3) instead of equa¬ 
tion (1). 

The right hand side term in equation (3) is shown in 
figure 5 to be small and to vanish after disorder averaging. 


Derivation of equation (7) 

We have from equation (5) 

c,(a:) = (fc;(0)e-^*r(14) 
Expanding ky{x) ~ fc„(0) -f xk'^{0) for |x| ^ I, yields 
Cpix) ~ 

~ (e-2“fc-'(o))(/fc(^(0)e-“'<(°)) (15) 

since fc(,(0) is not correlated to any function of fc„(0) 
because (fc(,(0)A:”(0)) = 0 for any positive integer n, 
when assuming random Gaussian impurities. The last 


term can now be expressed in terms of the distribu¬ 
tion function of the disorder potential ky, i.e., P(2x) = 
f dkyP{ky)e~'^^^'^'’ = Similarly we have 

iPp^x"^) = J dfc(,Pp(fc(,)A;(,e““^^», where Pp is the distri¬ 
bution function of the disorder potential ky. Hence we 
obtain equation (7). 

Disorder strength dependence around the resonance 

We used two different techniques to obtain a random 
potential characterized by a correlation length 1. In fig¬ 
ures 1-3 we used the sum of Gaussian impurities with 
random amplitudes and located at random sites. The 
number of impurities scales as 1/1 in order for the stan¬ 
dard deviation of the potential to be independent of 1. 
The advantage of this potential is that computing the bi¬ 
nary correlator is very simple as seen in equation (2). In 
figures 4 and 5 we used a random potential obtained by 
smoothing an uncorrelated potential over I neighbors us¬ 
ing the local regression smoothing process with tri-cube 
weight functions [35]. This produces a smooth poten¬ 
tial with a characteristic correlation length 1. Here we 
have no simple expression for the binary correlator or 
its Fourier transform, which has to be computed numer¬ 
ically. However, the Fourier transform of this correla¬ 
tor has resonances where the Fourier transform vanishes, 
which corresponds to the resonances seen in figure 4. A 
typical realization is shown in figure 5 as well as the de¬ 
pendence of A on disorder strength. 




FIG. 5. Left: the blue line shows the spatial dependence of 
ky{x) for a particular disorder realization. The green line is 
the difference between the quantum and the classical momen¬ 
tum {x)[—idx — p\'il!{x)), while the red line is the dis¬ 
order averaged difference (1000 realizations) (x)[—idx — 
p]%l}{x)). Here we used fco = 6.3, I = 1.15 and at = 0.9. 
Right: The Lyapounov exponent as a function of the disorder 
strength ak around a resonance {I = 0.85). The blueish lines 
are the numerical A from the transmission, the blueish sym- 
blols are from our approach (equation (6)), while the red lines 
are from the perturbative approach (2Au,). For I — 0.85 the 
perturbative approach gives Am = 0 regardless of disorder. 









